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We study orthogonal polynomials with weight exp[— NV(x)], where V(x) = 5^fe=i a 2kX 2k /2k is a 
polynomial of order 2d. We derive the generalised Freud's equations for d — 3, 4 and 5 and using 
this obtain = h^/h^-x, where h^ is the normalization constant for the corresponding orthogonal 
polynomials. Moments of the density functions, expressed in terms of 7? M , are obtained using Freud's 
equation and using this, explicit results of level densities as N — ► oo are derived. 
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Universality in random matrix theory 0, [Tl| - [l3| has led people to study orthogonal 0, 0, [j| and skew-orthogonal 
polynomials [3( in great details. However, in the process, the non-universal level densities are neglected inspite of the 
Q4 possibility of its direct application in various physical systems. In this context, we study level densities of a class of 
(Jh \ non-Gaussian random matrix ensembles and thereby develop the theory of orthogonal polynomials. 
» ■ Orthogonal polynomials are defined as 
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PACS numbers: 02.30.Gp, 05.45.Mt 



INTRODUCTION 



P n (x)P m (x)w(x)dx = h n 5 nm , n,meN. (1.1) 

We study orthogonal polynomials with weight function w(x) — exp(— NV(x)), where 

d 

Q2kX /(2k), Q2d>0- (1.2) 



fc=l 



Here, we make a numerical analysis of orthogonal polynomials corresponding to d — 3, 4 and 5. We derive the 
corresponding Freud's equation and calculate i? M = h^/h^-i. We observe interesting patterns in the behavior of 

Once we have an understanding of we use these results to obtain level densities of non-Gaussian ensembles of 
random matrices. We know that variation of the first n-eigenvalues of a random matrix can be studied by the n-point 
correlation function, R^f\xi, ...,x n ) which is defined by 

>< 



R^(xi,...,x n )= N ' [ dx n+ i...dx N Pp iN (xi,...,x N ), (1.3) 



where /3 — 1,2,4 correspond to ensembles of random matrices invariant under orthogonal, unitary and symplectic 
transformations. This allows us to find the probability density of the n eigenvalues at x±,...,x n , irrespective of the 
eigenvalues at x n+ i...xpj. R\^'(x) is the level density, which, for f3 = 2 can be written as [l9l - l22j 

JV-l 

Rf\x) = Y J {h^[P,{x)fe- Nv ^. (1.4) 

To calculate R\ (x) as N —t 00, the standard method is to use the Christoffel Darboux formula and the asymptotic 
results of orthogonal polynomials. The latter is not always available for general polynomial potential inspite of some 
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serious contributions from several authors [2J, Il5l4l8| on the asymptotics of orthogonal polynomials with V(x) = x 2d 
using the Riemann Hilbert technique. In this paper, we use the method of resolvent to obtain the level densities as 
TV —5- oo.This needs an understanding of moments Mfc defined as 

M k = [ x k Rf\x)dx, fceN. (1.5) 

Jr 

This is derived using the values of i? M using generalised Freud's equation, which we derive independently. Using 
this, we obtain the corresponding level densities. This gives us a good understanding of the origin of multiple band 
formation in the level densities in polynomial potential. 

The paper is organized as follows: In section 2, we study the d — 3 case and observe the behaviour of for 
different values of a^. Section 3 and 4 deal with d = 4 and d = 5 results. This is followed by our concluding remarks. 



2. d = 3CASE 



A. Freud's equation 

Orthogonal monic polynomials with even weight satisfy a recursion relation [l[ 

xP» = P»+i + Vh. (2.1) 

where i? M = hn/h^—i, for (i > 1 and i?o = 0. 

A major development in the study of quartic weight (d = 2 in eq. Ql.2p ) polynomials (25l - [27| was the following 
recursive equation in i? M due to Q. 

H + 1 = NR^ +1 [a 4 (R^ +2 + R„ +1 + R„) + a 2 \. (2.2) 
Now, we derive a similar Freud's equation for sextic potential, i.e. d = 3 in eq. dl.2p . We use the identity 

dx[P^ +1 (x)P^x)e~ NV ^Y = 0. (2.3) 
Using Pu{%) = + . . . and the orthonormality condition p. II) . we get 

~ NV(x) ][PUi( x ) P ^ x ) + P„+i(x)P^dx + J Nlaex 5 + a 4 x 3 + a 2 x]P t2+1 (x)P^x)e' NV{x) }dx = 
This gives us 

= J Na 6 x 5 P fl+1 (x)P^x)e- NV( ^dx}+ J Na 4 x 3 P^ +1 {x)P^x)e- NV ^dx}+ J Na 2 xP^ +1 {x)P^{x)e- NV ^dx} 
Using (|2.ip . we obtain 

(i+l = NR ll+ i[a & (R fi+2 (R tl + Rfj,+i + R^+2 + Rn+3) 

(2.4) 



Rfi+iiRfi + Rp+i + Rfi+2) 



+ Rij,(Rfj,-i + Rfi + Rfj,+i)) 
+ 04(^-1 + R» + R^+i) + 09]. 

Here we note that the corresponding Freud's equation is cubic in nature thereby giving rise to oscillatory behavior. 



B. The Kfi plot 

For the d — 2 case, two main features were observed in the R^ plot from the original Freud's equation: A two band 
structure formed by an oscillation between two values, converging to a single band. 

In the sextic case, the two band and single band structures reappear, however a new, more chaotic structure is also 
seen, appearing either between one band and one band or one band and two band structures. Henceforth, it is termed 
as a "transient structure". 
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FIG. 1: Rft plot modulo 3 for d = 3,a 6 = 1, 04 = -2.5, a 2 = 1, JV = 500 using eq. JO} 

1. Single band structure 

For the single band structure, all the terms in the Freud's equation are equal to each other. Thus, we obtain a A 
dependent cubic equation (where A = fi/N) , solving which we obtain one real solution 



D± = y [{2a\ - 10a 6 a 4 - 100a§A) ± J {2a\ - 10a 6 a 4 - 100a§A) 2 - A{a 2 6 - H^)3 



(2.5) 



04 



10a6 30a6 30ag 



(2.. 



which gives the value of i? M for A values where single band exists. 



2. Two band structure 

Solving the Freud's equation assuming that two bands are formed (as seen), i.e. 
Ao = Ro = i?2 = ^4 = ••• 
A\ = Ri = R3 = R5 = ... 
for N ^> n, we get 

Aa+Ai = ~°4 + y/aj - 4a 2fl6 
2a 6 

It has been numerically verihed that the bottom band (Al) tends to 0, and we find that A\ oc l/N and A% oc l/(a4) 2 



3. Transient structure 



When the transient structure is divided modulo 3 into three bands (&o, &i and 62 in fig. [p, it is seen that each of 
the three separate bands continuously oscillate, and converge to a common value. 
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The sum mod 3 for the duration of the transient structure oscillates above the value 



bo + bi + b 2 



— 04 + \Ja\ - 4a 2 a 6 
2a 6 



(2.8) 



where b , bi and b 2 correspond to three consecutive values from each of the bands. 



C. Critical a 4 's 

On analyzing the roots of the sextic potential, we obtain the critical value of 04, denoted here by a^ c , where the 
structure of the plot changes. We find the points at which potential plot touches 0, and obtain 



«4r 



48 



(2.9) 



In the i? M plot, when (i) 04 < a^ c , we observe a two band structure, followed by a transient structure, which 
converges to single band. (ii)<Z4 > a4 C , we see a one band structure, followed by a transient structure, which converges 
to single band. 

We now analyze the behavior of as 04 approaches a^ c (fig. [2]). It is observed the transient structure resolves 
into three distinct bands near a^ c - Exactly at the critical value, two of these bands coincide to form an upper band, 
and the third band forms a lower band, creating a pseudo two band structure. 
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FIG. 2: plot as a 4 approaches a 4c for ae=l, 02=1, N=500 
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D. Level Density 

In this subsection, we derive the level density using the method of resolvent 0, [Io[ as N — > oo. The result is 
expressed in terms of moments Mj, (|1.5I) which are derived using results from the Freud's equation. We then compare 
the result with the level density for TV = 30 using (|1.4[1 . 




i 



FIG. 3: Level density plots for d = 3, a§ = 1, <i4 = —3, 02 = 1 iV = 30. For the theoretical plot, calculated moments are 
M 2 = 62.536, M 4 = 164.770. 

Here, we recall Qthat we use the scaling V(x) — > V(x)/2 to obtain the corresponding results. 



V'(z)-V'(x) p(2) 



z — X 



R\'{x)dx~N z 



V'(x) 



N [a 6 (x 4 N + x 2 M 2 + M 4 ) + a 4 {x 2 N + M 2 ) + a 2 N] - N 2 x 2 



2 2 ( CLQX + d 4 X + 0,2 



Finally, 



TrR (2) (x) 



N 



a 6 M 4 + a 4 M 2 
N 



\ o \( 2 ^M 2 \ 1 , 4 2 S 2 

+ a 2 I + a; la 6 a; H — V a± J — - [clqx + a 4 x + a 2 ) 



(2.10) 



( 2) 

Now that we have derived R\ (x) in terms of M 2 and M4, we would be interested to calculate them using Freud's 
equation. We use 



x k P 2 (x) 

w(x)dx 



= EE 

= E c - 



w[x)dx 



(2.11) 



where are coefficients which can be expressed in terms of R^. A few typical examples are 



N-l 

M 2 = J2 (^M+l + R ») > 
M=0 



(2.12) 
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and 



N 

M4 = (i?^ + R^ + i + 2R[ i R fJ _+i + + R^R^-i) , 



N 

Me = (Rfj,~2Rfi-iRft + Rft-iRfi + ZRfi-iRjt + R^ + ^R^-iR^R^+i 



(2.13) 



+3Rj l R fl+ i + 3R fl R^ l+1 + + 2R ll R fJi+ iR tl+ 2 + 2Rj l+1 R fJ/+2 
+R fJ , +1 Rf l+2 + R^+xR^Rfi+a)- 



(2.14) 



The expression for Ms and higher moments are extremely cumbersome but can be easily calculated using the 
aforementioned algorithm. 



3. d = 4CASE 
A. Freud's equation 
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FIG. 4: d = 4 case plot for as = I, ag = —5, 04 = 6, ei2 = — 1, iV = 200 
As in the d = 3 case, we start with the identity 

dx[P^ +1 (x)P^x)e- NV ^]' = 0, 



(3.1) 
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where V(x) is defined as in eg. (|1.2|l . but with d = 4. Using the recursion relation for orthogonal polynomials, we get 

Ai+4 fi+3 M-2 

A 1 + 1 = NR^+% [a 8 (i? AI+2 i? AI +3 ^ -R; + ^ + -Ri 

i— i—^L i—fi 

+ Rfi+2Rfi ^ ^ + Rn+iRn+2 / ] Rj + -Ru+1 -R< + R^Rn+i ^ 

i=fi— 1 i— fi i=fi i=fj,—l 

)i+2 n+1 n+1 

+ RfiRfi+i^j R4 + RfiRn-i ^ Ri + R 2 ^ 2_j (3*2) 

i= i— fi — 2 i— fi — 1 

M+3 (U+2 

+ Q6(-R/j+2 -R» + -R/j+i Rj + Rp 

|U+1 

+ a 4 ^ i?j + a 2 ] 

Here we note that due to the non-linear nature of the Freud's equation, we observe oscillations in the solution for 
R^. These oscillations can be seen when the plot is divided modulo 4 into 4 residual bands (60,61,62 and 63 in fig. @| 



B. Level density 



Using the obtained moments and the formulation for finding level density (sec. 12 Dl) . we derive the function for 
R { i\x) for the d = 4 case (13. 3[) . The plot of iZ^(ac) obtained using this function is compared with the R^\x) 
calculated from (11.41) in fig. [5] 

We note that irregularities in the form of small oscillations around the expected value are seen near the peaks. This 
is because the numerically calculated level density is for a finite value of N. These oscillations gradually disappear as 
the value of N increases. 
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FIG. 5: Level density plots for d = 4, as = 1, ae = —5, 04 = 6, a? = — 1, JV = 20. For the theoretical plot, calculated moments 
are M 2 = 43.475, M 4 = 134.555, M 6 = 438.400. 
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4. d = 5 CASE 
A. Freud's equation 

As in the d = 3 case, we start with the identity 

dx[P^ +l {x)P^x)e- NV ^\ = 0, (4.1) 

where V(x) is defined as in eq. (jl.2l) . but with d — 5. From here, one can understand that finding the Freud's equation 
is algorithmic in nature and we leave it to the reader to derive it explicitly. Here, we will show the generic plot of the 
i? M function. 
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FIG. 6: d = 5 case plot modulo 5 for oio = 10, as = -80, a 6 = 210, a 4 = -200, a 2 = 48, N = 50 



B. Level density 



Having derived the moments M2, M4, M6 and Ms using the general formulation given in sec. 12 D\ we use the 
derivation provided to obtain the function for R[ (x) for the d — 5 case (j4~2l> . The plot of R\ >{x) obtained using 

( 2) _____ 

this function is compared with the R{ (x) calculated from (|1 .4|) in fig. [71 

Once again, we note that small oscillations around the expected value are seen. This is because we are calculating 

(2) 

R\ (x) for a finite value of N, and these become smooth as N — > 00 
-1 2 
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FIG. 7: Level density plot for for d — 5, aio = 1, ag = —8, as = 21, 04 = —20, 02 = 4.8, AT = 20. For the theoretical plot, 
calculated moments are M 2 = 43.960, K'h = 136.008, M 6 = 460.375, Ms = 1625.995. 



5. CONCLUSION 

In this paper, we obtain the Freud's equation for polynomials with weight function exp[— NV(x)], where V(x) = 
Sfe=l o,2kX 2k /2k is a polynomial of order 2d. We derive the generalised Freud's equations for d = 3, 4 and 5. We 
observe limit cycle behavior of R^. We use these results and the method of resolvent to obtain the level densities of 
the corresponding random matrix models. However, this involves an explicit calculation of the higher moments which 
we calculate numerically and insert in the analytic results of the level density. It would be nice to obtain explicit 
results of these moments as was done for the quartic case. But we have failed in this investigation due to the complex 
nature of the Freud's equation. 

Here, we might recall that for d — 2, the Freud's equation is quadratic while for higher d, it becomes cubic 
(d = 3), quartic (d = 4) and so on. This results in oscillations in the function and hence studying the limit 
cycle behavior becomes increasingly complicated. Further investigation is needed to study these non-linear Freud's 
equations, specially in the context of integrability and hence the existence of Lax pairs. We wish to come back to 
these questions in a later publication. 
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